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Abstract We have studied the acoustic Faraday effect in superfluid 3 He up to sig- 
nificantly larger magnetic fields than in previous experiments achieving rotations 
of the polarization of transverse sound as large as 1710°. We report nonlinear field 
effects, and use the linear results to determine the Zeeman splitting of the imagi- 
nary squashing mode (ISQ) frequency in 3 He-B. 
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1 Introduction 

Superfluid 3 He is a p-wave, spin-triplet superfluid of great interest due to its un- 
conventional pairing symmetry. While numerous experimental probes exist, trans- 
verse zero sound has recently become an excellent tool to explore the order pa- 
rameter structure, after its prediction by Landau 1 and Moores and Sauls 2 , and 
discovery by Lee et al. 3 The coupling between transverse sound and the Imagi- 
nary Squashing Mode (ISQ), an order parameter collective mode, allows for spec- 
troscopy of the mode and its behavior in a magnetic field. The ISQ has total an- 
gular momentum 7 = 2 with five Zeeman sub-states. Transverse sound couples 
to the rtij = ±1 states, causing a rotation of the sound polarization in a magnetic 
field, called the Acoustic Faraday Effect (AFE). It is this effect that both proves 
the existence of transverse zero sound in the superfluid 3 and provides a sensitive 
probe into the magnetic field dependence of the ISQ.— 

Transverse sound interacts with the ISQ for acoustic frequencies near the mode 
according to the dispersion relation,— 
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Fig. 1 (Color online) ISQ (blue line) and pair-breaking (green dashed line) frequencies plotted 
against reduced temperature. A typical pressure sweep alters the sound frequency normalized 
to the zero temperature gap to follow the red line, and transverse sound only propagates in the 
grey shaded area. Inset: An acoustic trace versus pressure at H = 0.04 T, illustrating the cavity 
oscillations from changes in the sound velocity and the envelope from the AFE. The black arrows 
indicate AFE rotations of 90°, 270°, and 450°, from right to left. 

where CO is the sound frequency, q the wavevector, Vf the Fermi velocity, and Q the 
ISQ frequency, with a zero-field value 5,4 of ~ */ 12/5 times the weak-coupling- 

plus gap. - Aq — yj- ( 1 — X) (1 + -y-) /(l + X -j-) is the quasiparticle restoring force, 

1.F S F s o F s 

and A 2 - = -jj-X(l + -j-) /(l +X-^-) is the superfluid coupling strength, where 
the F£ are Landau parameters, and X is the Tsuneto function. 7 We model the field 
dependence of the dispersion relation over a wide range of frequency by simply 
replacing the denominator on the right-hand side of Eq.[T]with 

D 2 = CO 2 - X2 2 - 2 -q 2 v 2 F - mjAy eff H - B^ S H 2 - mjC^H 3 . (2) 

Here Qq is the zero-field ISQ frequency, y e ff ls the effective gyromagnetic ratio 
of 3 He, H is the external magnetic field, and the terms containing A, B, and C de- 
scribe the linear, quadratic, and cubic magnetic field dependence of the dispersion, 
respectively. The m j dependence of each term is simplified to reflect the fact that 
nij = ± 1 , making the linear and cubic terms either positive or negative but having 
no effect on the quadratic term. 



2 Experiments 

In these experiments we use a dilution refrigerator and demagnetization stage 
setup described previously 4 to reach ~ 600 /iK with a 3 He cavity formed by a 
transducer and a reflecting plate. To probe the Faraday rotation, we lower the 
pressure of the helium from ~ 6 to 3 bar, changing the frequencies of the ISQ 
and pair-breaking, which transform relative to the acoustic frequency of 88 MHz 
as shown by the trajectory of the red arrow in Fig. [U taking into account heating 
inherent in the experiment which raises the temperature from ~ 0.5 — > 0.6T C . This 
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Fig. 2 (Color online) Data taken from acoustic traces at all experimental fields, plotted against 
the shift: (a) Sound polarization rotation angle data, 6. Open circles are taken from the envelope 
amplitude, while solid circles correspond to minima in the acoustic traces, (b) Sound velocity 
data, c,. The data points at each shift mostly overlap. 



causes changes in D 2 and thus the sound velocity, c t , leading to oscillations in the 
transducer response, or cavity oscillations, such as those seen in the inset to Fig. 
Q] With a magnetic field present, the polarization rotation is seen as the overall 
envelope in the inset with minima indicated by the black arrows. These effects are 
governed by 



Vz « cos 9 sin 



/2dC0\ 



(3) 



where Vz is the transducer voltage, 9 is the AFE rotation angle, and d = 31.6 ± 
0. 1 fim is the cavity spacing.— 

Using Eq. [3] we can get both 9 and c t from our acoustic response data. By 
measuring the envelope we calculate 9 based on its amplitude, using the nodes as 

fixed angles 9 = n x 90°; n = 1,3,5 We also convert the period of the cavity 

oscillations to changes in the sound velocity using 



tt) 

1 Period = 2d— 
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After calculating an initial value of c t near the ISQ, we use the sound velocity 
changes calculated from Eq.|4]to extend that initial sound velocity over the entire 
data range. 
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Fig. 3 (Color online) Values of X^ ec [ calculated using Eq.[7]for fitting 6 to the dispersion using 
A and C as free parameters, with B as the upper bound from fitting c,. 



3 Results and Discussion 

We show our data in Fig. [2] plotted against the normalized difference in the square 
of the frequencies, (ft) 2 — Qq ) /ft) 2 . This scaling, which we refer to as the frequency 
shift, or just the shift, better reflects the changes in the dispersion caused by the 
changes in temperature and pressure during an experimental run. 9 shows strong 
field dependence, but c t appears to be nearly field independent at these fields. 
We relate this data to the dispersion by taking Eq. [T| with the right-hand side 
denominator given by Eq.|2] and using the following definitions:^ 



9 = 2d 5q, 



(5) 



C, =2(0/(q++q-), (6) 

where 8q= \q + — q- 1/2, and q± is obtained by solving Eq.[T]for q, setting mj — 
±1. Using Eqs. [5] and [6] we find that c t depends most strongly on the quadratic 
field term, while 9 depends on the linear and cubic terms, due to the cancelling 
effects of mj = ± 1 . 

As mentioned above, c, shows very little field dependence. In order to quantify 
this we group the data in frequency shift bins spaced by (ft) 2 — Qq) /ft) 2 = 0.005, 
and fit each bin to Eq. Q] as a function of field, using Eq. [2] with B as the only 
free parameter and setting A = C = 0. All fits are done using the Levenberg- 
Marquardt algorithm^i for nonlinear least-squares minimization. We evaluate 
these fits quantitatively by calculating X^~a, a measure of how well a model fits 
data with a normal sample distribution, given by 



^red 



(yi-0 



(7) 



where V = N — M is the number of degrees of freedom, N is the number of data 
points, M is the number of fit parameters, y, is the measured value, fi is the cal- 
culated value, and a, is the standard deviation of the measurement. For c t we 
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Fig. 4 (Color online) A, B, and C parameters for Eq.[2] made dimensionless by normalizing to 
appropriate orders of CO. (a) The linear field dependence A/ CO, red circles, (b) An upper bound 
on the quadratic field dependence B, green triangles, (c) The cubic field dependence, C ■ co, solid 
orange diamonds. The upturn near co = C1q is in a region of greater uncertainty due to the lack 
of the highest field data, as described in the text. 



calculate a, directly from the data, assuming zero field dependence, resulting in 
values of a ~ 2 m/s, while for 9 we estimate a = 10° from our measurements. 
For a fit either a or 2a from our data, X^ e A ~ 1 -25 or 5.5, respectively. 

The results of our c t fits only give us an upper bound on its quadratic field 
dependence, which constrains the B term such that it has little to no effect on 9. 
Consequently, to account for the non-linear dependence of 9 when we fit 9 to 
Eq. Q] as above, we must include the cubic C term in Eq. [2] as well as the A term 
as free parameters, using the B from c t as an input, giving an accurate fit to the 
data resulting in the X^^a values shown in Fig. [3j At the lowest and highest shifts 

Xf e A increases, which is reflected in larger error bars for A and C at these shifts. 
We attribute these larger error bars to a larger uncertainty in 9 due to a greater 
difficulty in precisely identifying the rotation angle in these regions. Our A, B, and 
C results are shown in Fig. [4] given in a dimensionless form by normalizing to 
appropriate orders of tt). 

The appearance of an upturn at low shifts in C can be attributed to the lack of 
high field 9 data at those shifts, due to the rotation angle changing too fast to be 
determined. At higher shifts, C is a fairly constant value. Our data thus shows a 
decreasing linear term, an increasing quadratic bound, and a constant cubic term 
with increasing shift. 

Although there are predictions 1 1 for the quadratic field dependence of the ISQ, 
our B values cannot be directly compared with them due to the insensitivity of c t to 
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Fig. 5 (Color online) Lande g-factor values as a function of pressure, for our data (red circle, 
error bars inside data point) and reanalyzed data from Davis et al. (blue diamonds).-!— 

the magnetic field. As our data provides only an upper bound on the magnitude of 
the quadratic term in the dispersion we can only say that previous measurements-^ 
seem to fall within that bound at the mode. There are currently no theoretical 
predictions for the cubic field dependence of the ISQ. 

By extrapolating our A data to ft) = £2q, we can apply the theory of Sauls and 
Serene 13 to extract the Lande g-factor of the ISQ. The g-factor is a measure of 
the magnitude of the linear Zeeman splitting of the ISQ, and is predicted— to 
depend on the magnitude of /-wave pairing interactions in 3 He. As the theory 
applies only in the region around ft) ~ Q, we cannot extract more than a single 
data point, which, using the relationship Ao = 2gft), where Ao is the extrapolated 
value of A to ft) = £2q, is shown as a red circle in Fig.|5j We also present reanalyzed 
data from Davis et al. as blue diamonds. 4 The reanalysis, which focuses on their 
extrapolation of their data to T = rather than to ft) = Qq, as well as a discussion 
of the applicability of the theory, will be discussed in a forthcoming publication.— 



4 Conclusion 

We have measured the splitting of the ISQ using transverse sound, finding that 
nonlinear field effects play a significant role at fields up to H = 0.11 T. Using a 
simple model for the effect of splitting on the dispersion we have quantified the 
field dependence up to cubic order. We have set an upper bound on the quadratic 
splitting, which is small at these fields, and have determined the linear and cubic 
splitting parameters, as well as the Lande g-factor of the ISQ. 
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